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RESUME 

La  mesure  de  1' angle  d'arrivee  en  un  point  a  l'interieur  d'un 
faisceau  laser  se  propageant  dans  la  turbulence  est  obtenue  par  la  demo¬ 
dulation  en  frequence  du  signal  heterodyne  resultant  de  1 ' interference 
d'un  faisceau  sonde  avec  un  faisceau  de  reference.  On  utilise  ensuite 
ces  donnees  pour  demontrer  que  1 'angle  d'arrivee  et  l'intensite  lumineuse 
sont  faiblement  correles  et  que  1 'angle  d'arrivee  obeit  a  une  distribu¬ 
tion  de  probability  normale.  Ceci  constitue  une  verification  directe 
de  deux  hypotheses  fondamentales  du  modele  de  propagation  dans  la  turbu¬ 
lence  mis  au  point  au  Centre  de  recherches  pour  la  defense,  Valcartier 
(CRDV) .  (NC) 


ABSTRACT 

V 

Measurement  of  the  local  angle  of  arrival  in  a  laser  beam 
traveling  in  turbulence  is  achieved  by  frequency  demodulating  the  hete¬ 
rodyne  signal  derived  from  the  interference  of  a  probe  and  a  reference 
beam.  The  resulting  data  are  then  used  to  demonstrate  that  the  statis¬ 
tical  correlation  between  the  angle  of  arrival  and  the  irradiance  is 
weak  and  that  the  random  angle  of  arrival  is  normally  distributed.  These 
results  constitute  a  direct  confirmation  of  two  fundamental  assumptions 
of  the  theoretical  model  of  propagation  in  turbulence  previously  devel¬ 
oped  at  the  Defence  Research  Establishment  Valcartier  (DREV) .\(U) 
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1.0  INTRODUCTION 

The  phenomenon  of  turbulence  has  an  important  impact  on  optical 
and  infrared  laser  applications  in  the  atmosphere.  It  affects  both 
high  and  low  power  beams  and  is  generally  more  detrimental  under 
clear,  sunny  conditions.  The  main  resulting  optical  effects  on  laser 
beams  are  scintillation  and  beam  spreading.  Scintillation  is  a  random 
modulation  of  the  received  irradiance.  The  amplitude  of  the  modulation 
rapidly  grows  to  a  rms  level  of,  the  order  of  or  even  greater  than, the 
average  irradiance  and  the  structure  of  the  fluctuations  is  character¬ 
ized  by  short  and  intense  irradiance  peaks,  typically  greater  than  5  to 
10  times  the  average  irradiance.  These  modulations  can  obviously  cause 
problems  in  communication  coding  or  lower  the  breakdown  thresholds  at 
high  power  densities.  Beam  spreading  is  an  average  phenomenon  whereby 
the  many  turbulent  eddies  along  the  propagation  path  act  as  scatterers 
that  laterally  diffuse  the  beam  radiation.  Although  the  average  scatter¬ 
ing  angle  remains  relatively  small,  there  results  an  important  reduction 
in  the  beam  power  density  at  propagation  distances  of  a  few  kilometers. 
Turbulence  can  indeed  transform  into  a  net  loss  the  gain  normally 
predicted  by  focusing  the  beam  on  the  receiver. 

We  have  developed  at  DREV  a  mathematical  model  that  predicts  the 
average  irradiance  and  the  irradiance  variance  profiles  of  laser  beams 
propagating  in  turbulence  (Refs. 1-3).  The  model  can  also  simulate  the 
important  practical  situation  where  adaptive  corrections  are  applied 
to  the  transmitted  beam  (Refs. 4-5).  The  principal  advantages  of  our 
approach  over  existing  theories  are  that  it  is  uniformly  applicable  at 
arbitrary  scintillation  levels,  that  the  solutions  are  easily  calculated 
numerically,  and  that  extension  to  media  with  a  varying  and/or  a  non¬ 
linear  average  refractive  index  (e.g.  thermal  blooming  in  the  presence  of 
turbulence)  is  possible.  The  predictions  are  well  corroborated  by  data 
measured  in  the  atmosphere  and  in  laboratory-simulated  turbulence 
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(Refs. 1-3).  However,  one  basic  assumption,  although  proven  consistent, 
has  not  yet  been  verified  directly.  This  assumption  states  that  the 
instantaneous  irradiance  and  the  phase-front  angle  (or  the  angle  of 
arrival  of  the  optical  rays)  are  only  weakly  correlated.  In  this 
report,  we  describe  an  experiment  designed  to  measure  simultaneously  the 
local  random  irradiance  and  phase  of  a  laser  beam  traveling  in  turbu¬ 
lence  from  where  the  hypothesis  of  weak  correlation  can  be  tested.  The 
data  obtained  are  also  used  to  verify  directly  a  second  but  less  funda¬ 
mental  hypothesis  which  assumed  normal  statistics  to  relate  the  third- 
and  fourth-order  statistical  moments  of  the  angle  of  arrival  to  its 
second-order  moments . 

Section  2.0  reviews  the  experimental  approach  and  Section  3.0 
discusses  the  principal  results. 

This  work  was  performed  at  DREV  between  January  and  December, 

1980  under  PCN  33B07,  Atmospheric  Propagation  of  Laser  Beams. 


2.0  EXPERIMENTAL  APPROACH 


The  principal  measurements  presented  in  this  report  were  obtained 
by  a  heterodyne  detection  method.  This  section  describes  the  experi¬ 
mental  apparatus  and  briefly  reviews  the  theoretical  equations  pertinent 
to  the  present  application.  Also,  this  section  discusses  how  and  under 
what  conditions  the  heterodyne  signals  are  analyzed  to  obtain  the  tur¬ 
bulent  variations  of  the  angle  of  arrival. 
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2.1  Experimental  apparatus 

The  objective  of  this  work  is  to  measure  simultaneously  the  random 
fluctuations  of  the  irradiance  and  the  angle  of  arrival  of  a  laser  beam 
propagating  in  turbulence.  The  irradiance  is  obtained  by  ordinary  inco¬ 
herent  detection  but  the  angle  of  arrival  to  be  determined  necessitates 
a  coherent  heterodyne  detection  method.  The  latter  consists  in  mixing 
on  a  detector  the  probe  beam  that  propagates  in  the  turbulence  with  a 
reference  beam  that  travels  outside  of  the  turbulence  in  a  homogenous 
medium.  This  constitutes  an  interferometer  which  allows  the  measurements 
of  the  relative  phase  between  the  two  beams.  Since  the  reference  beam 
has  constant  parameters,  the  method  yields  the  state  of  the  instantaneous 
random  phase  of  the  beam  in  turbulence  from  which  the  angle  of  arrival 
is  derived. 


FIGURE  1  -  Schematic  diagram  of  the  experimental  apparatus.  BS:  beam 
splitter;  BE:  beam  expander;  M:  front -surface  mirror; 

NDF:  neutral-density  filter;  A/4 :  quarter-wave  retardation 
plate;  RP:  right-angle  prism;  CBS:  cube  beam  splitter; 

D:  light  detector. 
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The  experimental  apparatus  is  illustrated  in  Fig.  1.  Both  the 
probe  and  the  reference  beams  are  extracted  from  the  same  50-mW  He-Ne 
laser  source.  The  reference  beam  is  reflected  from  the  front  surface  of 
a  glass  window  positioned  immediately  at  the  exit  of  the  laser;  it  has 
a  2-mm  diameter  and  contains  about  4%  of  the  total  power.  The  reflection 
from  the  back  surface  is  blocked  up.  The  transmitted  part  is  expanded 
to  a  25-mm  diameter  and  constitutes  the  probe  beam  which  is  deflected 
by  mirrors  into  the  turbulent  medium. 

The  turbulent  medium  is  a  scaled  simulation  of  the  atmosphere 
according  to  a  technique  described  in  Refs. 6-7.  In  short,  the  refrac¬ 
tive-index  fluctuations  are  produced  by  creating  an  unstable  vertical 
temperature  gradient  in  a  tank  filled  with  water.  This  is  simply  done 
by  heating  the  water  at  the  bottom  of  the  tank  and  cooling  it  at  the 
top.  The  tank  is  1.5  m  long,  0.6  m  deep,  and  0.4  m  wide.  The  charac¬ 
teristic  Kolmogorov  inertial  subrange  is  reasonably  well  verified  and 

-4  -1/3 

the  index  structure  constant  is  typically  10  m  .  The  propaga¬ 
tion  axis  is  along  the  length  of  the  tank  and  the  distance  can  be 
increased  by  folding  the  beam  as  many  times  as  required. 

The  measurement  method  for  the  angle  of  arrival  is  based  on  the 
frequency  demodulation  of  the  heterodyne  signal  obtained  by  mixing  the 
probe  and  the  reference  beams.  But,  this  approach  is  possible  only  if 
the  modulation  caused  by  the  fluctuating  angle  of  arrival  is  carried  by 
a  higher  frequency  signal,  which  is  done  here  by  shifting  the  optical 
frequencies  of  the  probe  and  the  reference  beams  with  respect  to  one 
another.  Considering  all  the  operating  conditions,  we  find  that  a 
convenient  shift  should  be  of  a  few  kilohertz.  Reference  6  gives  a 
review  of  several  techniques  of  optical  frequency  shifting  available 
in  this  range,  but  the  one  best  suited  to  our  needs  is  the  use  of  an 
electro-optic  crystal.  As  applied  in  the  present  experiment  and  illus¬ 
trated  in  Fig.  1,  the  technique  consists  in  letting  the  reference  beam 
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propagate  through  a  30-mm-long,  1  x  1  mm  crystal  of  LiNbO^  excited  by  a 
transverse  electric  field  rotating  at  a  frequency  n/2ir.  The  latter  is 
simply  produced  by  applying  two  sinusoidal  voltages  of  equal  amplitude 
(nominally  200  V)  and  in  quadrature  to  the  lateral  faces  of  the  crystal. 
The  principle  of  operation,  outlined  in  Ref.  6,  is  that  the  rotating 
electric  field  inside  the  crystal  has  the  property  of  accelerating  (or 
decelerating)  the  incident  circularly  polarized  field  of  the  laser  light 
beam,  thus  creating  an  upshift  Q.  (or  downshift)  of  the  angular  optical 
frequency  of  the  outgoing  light.  The  method  is  equivalent  to  a  rotating 
half-wave  plate  and  the  efficiency  quoted  in  Ref.  6  is  55%,  which  is 
well  corroborated  by  our  measurements.  Two  X/4  -  phase  retardation 
plates  are  used  to  transform  the  linear  polarization  of  the  reference 
beam  into  circular  polarization  and  vice  versa.  The  first  lens  focuses 
the  beam  inside  the  crystal  for  better  efficiency  and  the  second  one 
recollimates  it.  Following  the  frequency  shifter,  a  set  of  right-angle 
prisms  is  arranged  to  adjust  the  propagation  path  of  the  reference  beam 
to  that  of  the  probe  beam.  This  is  necessary  since  the  longitudinal 
coherence  length  of  the  laser  is  only  a  few  tens  of  centimeters. 

Finally,  the  two  beams  are  brought  together  or  mixed  by  a  50-mm  cube 
beamsplitter  as  shown  in  Fig.  1;  the  interfering  beams  are  then  detected 
by  a  silicon  phototransistor  mounted  behind  a  0.5-mm-diameter  aperture. 

2.2  Heterodyne  Detection 

The  reference  and  the  probe  beams  are  linearly  polarized  in  the 
same  direction  and,  as  is  customary  and  justified  (Ref.  7,  p.  93,  eqs. 
6.1  and  6.2)  in  the  treatment  of  atmospheric  turbulence,  it  is  assumed 
that  polarization  is  not  affected  by  propagation.  A  scalar  represen¬ 
tation  is  therefore  possible  and  the  electric  field  E  in  the  respective 
beams  is  denoted  as  follows: 

E  =  A  expf-iw  t  +  ikn  z+ikn  4>  1 ,  [1] 

p  P  P  op  p  op  pJ  ’  1  1 
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Er  =  Ar  exp[- 


10)  t  + 
r 


ikn  z 
or  r 


ik  n  $  ] 
r  or  rJ 


[2] 


where  A  is  the  amplitude,  u)  is  the  optical  angular  frequency,  k  is  the 
free-space  optical  wave  number,  c  is  the  speed  of  light  in  free  space, 
nQ  is  the  unperturbed  index  of  refraction  of  the  medium,  z  is  the 
geometrical  distance  of  propagation,  <f>  is  the  phase  shift,  i.e.  the 
difference  between  the  optical  and  the  geometrical  thicknesses  of  the 
propagation  medium,  and  i  -yfr.  The  subscripts  r  and  p  refer  to  the 
reference  and  the  probe  beams  respectively.  The  two  beams  are  mixed 
and  aligned  parallel  to  each  other  on  a  photodetector  which  thus  gives 
a  current 


(E. 


E  )  (E  + 
r  P 


E  )’ 

r 


[3] 


where  the  superscript  *  denotes  a  complex  conjugate.  Substituting 
eqs.  1-2  for  E  and  E^  in  eq.  3  and  assuming  that  ft  =  wr  - 


smaller  than  either  w  or  u  ,  we  find 

r  p 


io  is  much 
P 


j  A2  +  A2  +  2A  A  cos[S2t  +  k(n  4>  -n  <j>  )+k(n  z  -n  z  )].  [4] 

J  r  p  r  p  op  p  or  r  op  p  or  r  J 


op  p 


For  detection  at  a  given  fixed  point,  the  parameters  A^,  <J>r ,  z^,  z  , 

n  ,  n  ,  SI  and  k  are  constant.  Hence,  eq.  4  can  be  more  conveniently 
or  op  n 

rewritten 


A2  + 
r 


A2 

P 


2A  A  cos  [Sit 
r  p 


+  kn  4>  +  cstl 

°P*P 


-  ...  - -  - f 


[5] 


UNCLASSIFIED 

7 


where  cst  is  a  constant  that  regroups  the  remaining  invariant  parameters, 
The  quantities  of  interest  to  the  present  s:udy  are  A2  and  <j>  ,  respecti¬ 
vely  the  irradiance  and  the  phase  of  the  probe  beam  that  has  traveled 
through  the  turbulence . 


The  requirement  of  the  following  analysis  is  to  derive  both 
quantities  A2  and  4^  from  the  single  heterodyne  signal  given  by  eq.  5. 

A  simple  examination  of  eq.  5  shows  that  by  choosing  the  frequency  shift 
12/2rr  much  greater  than  the  highest  frequencies  of  the  turbulent  irra¬ 
diance  A2,  it  is  possible  to  separate  electronically  the  first  two  terms 
from  the  last  one.  Indeed,  if  the  cutoff  frequencies  are  properly  set, 
a  low-pass  filter  will  yield  a  signal 


e.  «  A2  +  A2  , 

lp  r  p 


and  a  high-pass  filter,  a  signal 


e,  «  A  A  cosfflt  +  kn  4>  +  cst]. 

hp  r  p  L  op  p  1 


Since  A^  is  constant,  it  is  easy  to  subtract  its  contribution  from 
eq.  6  and  we  thus  obtain  a  signal  proportional  to  the  fluctuating 
irradiance  A2  of  the  probe  beam  at  the  position  of  the  detector. 


The  electronic  signal  given  by  eq .  7  is  a  signal  at  frequency 
fl/ 2n  modulated  in  amplitude  by  and  in  frequency  by  4  .  This  is  well 
illustrated  in  Fig.  2,  where  two  time  records  of  the  high-pass  filtered 
signal  are  reproduced:  a  first  one  at  a  slower  sweep  rate  to  show  the 


Cm  V  J 
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FIGURF.  2  -  F.xample  of  a  detected  heterodyne  signal  modulated  in  ampli¬ 
tude  and  frequency  by  the  turbulence;  propagation  distance 
z  =  1.62  m.  Trace  (al ;  slower  sweep  rate  to  exemplify  the 
amplitude  modulations.  Trace  (bl ;  faster  sweep  rate  to 
exemplify  the  frequency  modulations. 
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amplitude  modulations  and  a  second,  at  a  faster  sweep  rate  to  show  the 
frequency  modulations.  Since  the  frequency  domain  of  is  much  smaller 
than  ft/2ir,  the  characteristic  amplitude-modulation  periods  extend  over 
many  oscillations  at  frequency  ft/2rr.  Therefore,  it  is  possible,  in 
principle,  to  frequency  demodulate  the  signal  of  eq.  7  independently 
of  the  amplitude  modulations  and  thus  obtain  the  time  derivative  of 
the  phase  <(>  .  This  is  the  approach  chosen  in  the  present  experiment. 

The  apparatus  used  to  perform  the  frequency  demodulation  is  a 
Doppler  Signal  Processor  designed  for  laser  anemometry  and  built  by 
DISA  ELEKTRONIK  A/S.  Omitting  the  details  of  the  electronics  which  can 
be  found  in  Ref.  8,  this  apparatus  delivers  an  output  voltage  propor¬ 
tional  to  the  instantaneous  frequency  of  the  input  signal  it  tracks, 
i.e.  proportional  to  the  time  derivative  of  the  cosine  argument  of  eq.  7. 
The  output  voltage  of  the  Doppler  Signal  Processor  is  thus 

d  4> 

e,  =  Kffl  +  kn  — ]  ,  [8] 

dp  1  °p  dt 


where  K  is  the  calibration  constant.  This  output  voltage  e^  will, 
hereafter , often  be  referred  to  as  the  frequency-demodulated  heterodyne 
signal.  The  frequency  ft/ 2tt  is  adjustable  but,  for  most  results  reported 
here,  it  was  set  at  4.S  kHz.  The  Signal  Processor  was  operated  in  the 
frequency  range  of  2.25  to  15  kHz  which  gives  an  upper  tracking  fre¬ 
quency  limit  for  d^/dt  of  about  100  Hz  (Ref.  8).  This  limit  is  quite 
sufficient  for  the  type  of  turbulence  studied  here  as  the  results  will 
show.  Finally,  the  tracker  sensitivity  in  this  range  is  k  =0.100  ±0.002 
mV/rad/s  determined  from  the  slope  of  an  experimental  voltage-frequency 
calibration  curve. 
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2 . 3  Angle  of  Arrival 

The  quantity  derived  from  the  analysis  of  the  tracker  output  is 
d$  /dt .  However,  our  purpose  is  to  measure  the  instantaneous  angle  of 
arrival.  Following  our  propagation  model  developed  in  Refs.  1-3,  the 
angle  of  arrival  v,  at  point  (z,  _r)  is  defined  as  the  vector  angle 
subtended  by  the  geometrical  ray  which  passes  through  the  point  (.z,  _r)  ; 
z  is  the  propagation  distance  and  jr,  the  position  vector  in  the  plane 
perpendicular  to  the  propagation  axis.  From  Refs.  1-3,  we  thus  have  in 
the  notation  of  this  report 


v(=.  r)  »  7*p  (=,  r), 


[9] 


where  £  is  the  gradient  operator  with  respect  to  the  coordinates  of  _r. 
Hence,  we  need  to  relate  difi^/dt  to  .  The  sign  ~  is  used  in  eq.  9 
instead  of  =  because  the  definition  of  the  phase  <{>  differs  slightly 
from  that  of  Refs.  1-3.  On  the  one  hand,  the  phase  4>  used  in  Refs.  1-3 
satisfies  the  eikonal  equation  of  geometrical  optics,  i.e. 


—  +  -  (V^)2  =  (M-n  1  ,  110] 

3z  2 


where  N  is  the  instantaneous  index  of  refraction  of  the  medium  and 
where  the  paraxial  approximation,  valid  for  the  propagation  in  atmosphe¬ 
ric  turbulence,  has  been  made.  On  the  other  hand,  the  phase  $  of  eq.  1 
is  the  sum  of  $  plus  a  contribution  due  to  diffractional  effects  since 
the  amplitude  in  eq.  1  is  taken  as  real.  However,  diffraction  contributes 
only  slightly  to  the  fluctuations  of  $  and  from  here  on,  we  will  set 
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The  turbulent  angle  of  arrival  _v  =  is  generally  small,  typi¬ 
cally  smaller  than  1  mrad ,  even  in  very  strong  scintillation.  Hence, 
the  nonlinear  term  in  eq.  10  can  be  neglected  with  good  accuracy  and 
we  derive 


dip 

dt 


♦  W  — 
3z 


+ 


U-VN 


[11] 


where  W  and  ]J  are  respectively  the  fluid  velocity  components  along  and 
normal  to  the  propagation  axis.  The  first  term  on  the  right-hand  side 
of  eq.  11  is  negligible  since  it  represents  the  temporal  index  variations 
only,  i.e.  after  the  fluctuations  due  to  spatial  movements  of  the  fluid 
have  been  left  out.  These  temporal  variations  are  generally  much 
slower,  so  that  j  3N/3t  |  <<  | JJ*V,N  [  .  Similarly,  the  second  term  is  also 
neglected  since  the  rapid  fluctuations  of  3N/3z  are  averaged  out  by 
the  integration  over  z,  which  is  not  the  case  for  the  transverse  deri¬ 
vative  VN.  Moreover,  as  will  be  seen  in  the  following  discussion,  the 
velocity  component  W  in  the  regions  that  contribute  most  to  the  inte¬ 
gral  is  much  smaller  than  JJ.  Hence,  eq.  11  is  simplified  to 


chp 

dt 


U*Vn  dz 

iV  <-w/ 


[12] 


The  unstable  temperature  convection  that  characterizes  boundary 
layers  with  low  velocity,  as  is  the  case  in  the  atmosphere  and  in  our 
turbulence  simulator,  has  the  structure  of  vertical  jets  or  plumes. 

The  warmer  and  more  turbulent  fluid  rises  in  these  jets  while  the 
cooler  and  less  turbulent  fluid  sinks  in  more  uniform  and  less  localized 
currents.  This  phenomenon,  briefly  discussed  in  Refs. 9-10,  helped 
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explain  the  observed  intermittency  of  the  temperature  fluctuations 
recorded  in  our  simulator.  Three  important  properties  of  this  struc¬ 
ture  of  unstable  convection  are  used  to  further  simplify  eq.  12.  Firstly, 
the  level  of  temperature  or  refractive-index  turbulence  is  higher  inside 
the  jets.  Although  no  systematic  measurements  of  these  differences  have 
been  carried  out  in  our  simulator,  the  results  of  Refs. 9-10  indicate 
that  a  factor  of  3-S  is  probably  correct.  Secondly,  the  fluid  velocity 
inside  the  jets  is  greater  than  the  sinking  velocity  by  a  factor  of  4 
to  5  (Ref.  11,  p.  147).  Finally,  the  rising  velocity  is  approximately 
constant  over  most  of  the  jet  cross  section.  Therefore,  the  contribu¬ 
tions  to  the  integral  of  eq.  12  come  mainly  from  the  optical-path  seg¬ 
ments  that  pass  through  these  jets.  Hence,  if  we  assume  that  the  jets 
have  an  average  rising  velocity  U,  we  can  approximate  eq.  12  as  follows: 

z 

^  U  f  —  N  dz  ,  [13] 

dt  J  3y 

o 


where  y  is  the  vertical  coordinate.  In  other  words,  the  horizontal 
motion,  including  the  component  W,  is  less  important  and  contributes 
little  to  d<|>/dt.  Finally,  from  the  definition  of  the  angle  of  arrival 
v  =  £<t>  and  the  approximation  of  linearization,  we  obtain  from  eq.  19 

z 

X(z>  r)  =  j  VN  dz  .  [14] 

*  o 


Therefore,  comparing  eq.  13  to  eq.  14,  we  find 
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£4- 

dt 


[15] 


where  is  the  component  of  v_  in  the  vertical  direction. 

Equation  IS  relates  the  temporal  derivative  of  the  turbulent  phase 
fluctuations,  measured  by  demodulating  the  heterodyne  signal,  to  the 
vertical  component  of  the  angle  of  arrival.  This  equation  was  derived 
through  a  series  of  approximations  difficult  to  verify  individually. 
However,  the  data  analysis  will  show  that  eq.  15  is  consistent  and  can 
indeed  be  used  for  the  present  purpose  of  calculating  the  statistical 
correlation  between  the  irradiance  and  the  angle  of  arrival,  and  the 
probability  distribution  of  the  latter. 

2.4  Problem  of  Vibrations 


In  the  expression  for  the  heterodyne  signal  (eq.  5),  we  have 
collected  in  a  single  additive  constant  all  the  contributions  related 
to  the  geometrical  propagation  paths  z^  and  z^.  This  is  justified  in 
principle  since  the  detection  is  made  at  a  fixed  point.  In  practice, 
however,  there  are  oscillations  in  the  relative  path  difference  which 
can  obscure  the  useful  turbulent  phase  signal .  From  the  expected 
magnitude  of  the  turbulent  phase  fluctuations,  we  realize  that  oscilla¬ 
tions  in  (Zp  -  zr)  of  amplitude  as  small  as  1  ym  are  sufficient  to  do 
this.  Such  oscillations  are  difficult  to  eliminate  completely. 


To  minimize  this  problem,  all  the  components  of  the  experimental 
setup,  including  the  turbulence  simulator,  were  mounted  on  a  vibration- 
isolated  table.  The  improvement  over  the  nonisolated  conditions  is 
substantial  but  insufficient  to  consistently  bring  the  noise  level  due 
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to  the  vibrations  well  below  the  signal  level.  Figure  3  shows  an 
example  of  a  measured  spectral -density  function  of  the  frequency 
demodulated  heterodyne  signal.  This  density  function  was  calculated 
with  a  HONEYWELL  SAICOR-42A  Correlation  and  Probability  Analyzer  inter¬ 
connected  with  a  HONEYWELL  SAIC0R-470  Fourier  Transform  Analyzer;  the 
electronic  signal  was  sampled  at  a  rate  of  500  Hz,  and  2^  points  were 
used  for  the  computations.  We  observe  in  the  spectrum  of  Fig.  3  an 
important  peak  centered  at  about  32  Hz.  This  contribution,  which 
accounts  for  a  large  fraction  of  the  total  power,  is  caused  by  the 
vibrations.  Extensive  investigations  indicate  that  the  position  of  the 
peak  is  practically  independent  of  the  various  noise  sources  in  the 
laboratory,  pointing  to  a  phenomenon  of  resonance.  However,  this  contri¬ 
bution  cannot  be  removed  by  simple  electronic  filtering  since  it  occurs 
well  inside  the  expected  frequency  bandwidth  of  dij>/dt,  the  signal  of 
interest.  Moreover,  satisfactory  elimination  by  mechanical  isolation 
would  require  extreme  care  especially  because  of  the  large  number  of 
components  involved.  Therefore,  a  modified  approach  to  the  measurements 
was  designed  to  minimize  the  effects  of  vibrations  while  satisfying  the 
requirements  of  this  study. 

The  vibration  frequencies  observed  in  the  typical  spectral -density 
function  plotted  in  Fig.  3  are  sufficiently  low  to  suggest  that  the 
responsible  optical  components  oscillate  in  their  fundamental  mode. 

This  means  that  the  components  vibrate  as  a  solid  block  and  that  higher- 
order  rotation  or  torsion  modes  are  unlikely.  If  this  is  correct,  the 
vibrations  should  have  the  same  effect  irrespective  of  the  measurement 
position  across  the  beam.  Therefore,  if  the  heterodyne  signals  were 
recorded  simultaneously  at  two  points  on  the  beam  cross  section  and  if 
the  resulting  frequency-demodulated  signals  were  subtracted,  the  vibra¬ 
tion  contributions  should  cancel  out.  This  is  the  technique  we  have 
selected  and  Fig.  4  illustrates  how  it  is  implemented.  The  reference 
beam  is  split  into  two  parallel  beams  and  both  portions  are  mixed  with 


m 


Spectral  density  CdB/Hz) 
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Frequency  CHz) 


FIGURE  3  -  Example  of  the  measured  spectral  density  of  the  frequency- 
demodulated  heterodyne  signal  detected  at  one  point.  The 
peak  centered  at  about  32  Hz  is  due  to  the  effects  of 
mechanical  vibrations. 
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R  :  REFERENCE 
BEAM 


P  :  PROBE  BEAM 


FIGURE  4 


Schematic  diagram  illustrating  the 
heterodyne  detection  3t  two  points 
section.  BS :  beam  splittei;  CBS: 
Dl,  D2:  light  detectors. 


technique  for  simultaneous 
in  the  probe-beam  cross 
cube  bean,  splitter; 
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the  25-mm-diameter  probe  beam.  The  interfering  beams  are  then  detected 
simultaneously.  The  measurement  points  are  15  mm  apart,  positioned  at 
approximately  equal  distance  from  the  beam  axis. 

Ideally,  the  technique  would  require  two  Doppler  Signal  Processors 
to  demodulate  simultaneously  the  two  heterodyne  signals.  Since  only 
one  was  available,  the  analysis  was  performed  in  two  consecutive  opera¬ 
tions  as  shown  in  Fig.  5.  In  the  first  operation  (Fig.  5a),  one  of  the 
high-pass  filtered  heterodyne  signals  was  recorded  on  one  channel  of  a 
magnetic  tape  recorder  while  the  second  signal  was  demodulated  and  the 
output  of  the  processor  recorded  on  a  second  channel .  Logarithmic 
amplifiers  were  used  to  reduce  the  turbulent  amplitude  modulations. 

Also  recorded  in  this  first  operation  were  the  low-pass  filtered  signals, 
i.e.  the  irradiance  fluctuations  at  both  points.  In  the  second  opera¬ 
tion  (Fig.  5b),  the  tape  was  played  back  and  the  recorded  heterodyne 
signal  was  demodulated  in  its  turn.  The  output  was  then  subtracted 
from  the  previously  recorded  demodulated  signal.  The  resulting  voltage 

is  thus  proportional  to  the  required  function  (dij>  /dt-d<(>  /dt)  where 

1  2 

the  subscripts  1  and  2  refer  to  the  two  points  of  measurement. 

Figure  6  shows  the  spectral  density  function  of  (d<)>  /dt-d<}>  /dt) 

1  2 

for  the  same  time  segment  as  for  the  corresponding  spectral  density  of 

d<t> ^ / dt  plotted  in  Fig.  3.  The  peak  at  32  Hz  has  almost  completely 

disappeared  suggesting  that  the  vibration  contributions  have  cancelled 

out.  To  further  illustrate  how  the  technique  works,  we  have  also 

reproduced  in  Fig.  7  time  recordings  of  d<f>  /dt  and  (d$>  /dt -dd>  /dt)  for 

1  12 

nonturbulent  conditions.  In  the  latter  case,  (d$  /dt-d<j>  /dt)  should 

1  2 

be  zero  everywhere.  We  find  that  the  signal  level  has  indeed  dropped, 
by  a  factor  of  about  5  in  rms  values,  and  that  the  periodic  oscillations 
have  been  eliminated  by  subtraction.  Hence,  the  method  eliminates 
efficiently  the  effects  of  mechanical  vibrations  and  modifies  only 
slightly  the  analysis  of  the  data.  The  residual  noise  is  attributed  to 
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(a) 


Cb) 


FIGURE  5  -  Block  diagram  of  the  electronic  processing  of  the  heterodyne 
signals.  (a)  recording  operation;  (b)  playback  operation. 
D1 ,  D2:  light  detectors. 
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FIGURE  6  -  Example  of  the  measured  spectral  density  of  the  difference 

between  the  frequency-demodulated  heterodyne  signals  detected 
simultaneously  at  two  points.  The  contributions  due  to  the 
mechanical  vibrations  have  cancelled  out  almost  completely. 


Example  of  cancellation  of  the  etteets  of  mechanical  vibrations 
Trace  (a) :  time  recording  of  the  frequency-demodulated  hete¬ 
rodyne  signal  detected  at  one  point.  Trace  (b) :  time  record¬ 
ing  of  the  difference  between  the  frequency-demodulated 
heterodyne  signals  detected  simultaneously  at  two  points 
separated  by  15  mm  in  the  probe -beam  cross  section.  The 
recordings  were  made  under  conditions  of  vanishing  turbulence. 
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the  electronic  processing,  in  particular  to  the  Doppler  Signal  Processor. 
The  following  sections  will  discuss  how  the  useful  signal  and  the 
processing  noise  are  separated  for  each  application. 


3.0  RESULTS 

3.1  Variance  of  the  Angle  of  Arrival 

From  eqs.  8  and  15  and  the  assumption  that  the  contributions  due 
to  the  vibrations  cancel  out  exactly,  the  difference  between  the  output 
voltages  of  the  Doppler  Signal  Processor  for  the  two  measurement  points 
r  and  r  is  given  by 


e  -  e 
1  2 


Kkn  U  [VVCE, )  -  V. 


(r  )]. 

n 


f 

[16] 


The  variance  of  (e  -  e  )  is  therefore  proportional  to  the  structure 
1  2 

function  of  the  angle  of  arrival  for  a  separation  distance  A_  =  r  -  r 
which  was  chosen  along  the  horizontal  axis,  i.e.  A  =  Ai.  From  the 
theoretical  results  of  Ref.  3  on  the  covariance  of  the  angle  of  arrival, 
we  then  have 


<(e  -e  )2>  =  K2k2n2  U2  <fv  (r  )  -  v  fr  ))2? 

12  o  y  ~i  y^2 


=  6.90 


K2k2n2U2C2z 
o  n 

-7T7T— 


{l  -  [1+8.16  A2/£2]‘1/6J  , 


[17] 
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where  the  pointed  brackets  <...>  denote  ensemble  averaging  and  where 

C  is  the  refractive-index  structure  constant,  £  is  the  inner  scale 
n  o 

of  turbulence,  and  z  is  the  propagation  distance. 

Equation  17  is  a  variant  of  eq.  15  which  relates  the  measured 
time  derivative  of  the  turbulent  phase  to  the  angle  of  arrival.  As 
discussed  previously,  eq.  15  is  central  to  the  present  analysis  but  it 
is  based  on  a  series  of  approximations  difficult  to  verify  individually. 
However,  eq.  17  now  provides  a  means  of  determining  its  consistency. 

To  do  this,  the  variance  of  (e  -e  )  was  measured  at  several  propagation 

12 

distances  leaving  all  the  other  parameters  constant.  Equation  17  shows 
that  these  data,  when  plotted  against  z,  should  fall  on  a  straight  line 
through  the  origin,  which  is  reasonably  well  confirmed  by  Fig.  8. 

The  contributions  due  to  the  processing  noise  that  are  subtracted 

from  the  data  of  Fig.  8  were  calculated  by  recording  the  heterodyne 

signal  of  one  of  the  detectors  along  with  the  demodulated  signal  from 

the  same  detector.  On  playing  back  the  tape  and  performing  the  same 

analvsis  as  for  <(e  -e  )2> ,  we  thus  obtain  the  residual  noise  level. 

1  2 

This  procedure  was  repeated  for  all  the  data  points  of  Fig.  S  and  we 
found  that  the  noise  was  constant  to  within  r  S'1.,.  All  the  variances 
were  measured  with  the  HONEYWELL  SAIC0R-42A  Correlation  and  Probability 
Analyzer  from  21  samplings  of  the  instantaneous  electronic  signals 
digitized  at  a  rate  of  5iU’  Hz. 

The  processing  noise  has  its  origin  in  the  strong  turbulent 
amplitude  fluctuations  of  the  heterodyne  signal  characteri zed  by  intense 
peaks  followed  by  periods  of  vanishing  amplitude.  If,  during  these 
periods,  the  signal  drops  below  a  preset  threshold,  the  Doppler  Signal 
Processor  breaks  its  servo  loop  and  holds  its  output  until  the  signal 
reappears  at  a  frequency  within  certain  iimits  of  the  last  measured 


FIGURE  8  -  Variance  of  the  difference  between  the  frequency-demodulated 
heterodyne  signals  detected  simultaneously  at  two  points 
separated  by  15  mm  plotted  versus  the  propagat^on^d^stance  for 
constant  turbulence  parameters  C  =  1.05  x  10"  m  and 
i  =  2  mm.  The  contributions  due  to  the  processing  noise 
were  subtracted  from  these  data.  The  error  bars  represent 
the  standard  deviation  of  the  measured  noise  level  for  the 
six  data  points. 
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Doppler  frequency  whereupon  tracking  is  resumed,  usually  accompanied  by 
a  rapid  change  of  the  output  voltage.  On  the  other  hand ,  if  the  thres¬ 
hold  is  set  too  low,  the  processor  may  track  a  false  noise  signal  and 
thus  give  an  erroneous  output.  The  problem  is  most  strongly  felt  in  the 
analysis  of  the  recorded  heterodyne  signal.  Because  the  dynamic  range 
of  the  recorder  is  limited,  it  is  difficult  to  accommodate  the  large 
turbulent  amplitude  modulations  and  we  find  that  the  signal  can  fall 
below  the  tape  noise  level.  Since  the  tape  noise  is  in  the  same 
frequency  range  as  the  frequency  ;./ 2n  of  the  heterodyne  signal,  the 
instrument  tends  to  track  the  noise  during  these  periods.  In  any  case, 
whether  they  are  due  to  tracking  dropouts  or  to  tracking  errors,  we 
observe  frequent  short  peaks  in  the  processor  outputs  that  are  correlated 
with  the  periods  of  near-zero  amplitude  of  the  recorded  heterodyne 
signal.  Through  careful  adjustments  of  the  recorder  dynamic  range  and 
of  the  Doppler  Signal  Processor  threshold,  we  can  minimize  these  errors 
but  we  cannot  eliminate  them  completely.  However,  since  the  duration 
of  the  peaks  are  much  shorter  than  the  inverse  bandwidth  of  the  useful 
phase  signal,  they  are  uncorrelated  with  the  latter.  Hence,  their 
contribution  to  the  measured  variance  of  the  processor  output  can  be 
subtracted  with  good  accuracy  even  though  it  constitutes  a  relatively 
large  fraction  of  the  total  variance.  In  the  processing  of  the  direct 
or  nonrecorded  heterodyne  signal,  the  error  peaks  can  be  almost  comple¬ 
tely  eliminated  by  proper  threshold  adjustments.  Therefore,  the  residual 
noise  comes  mainly  from  channel  No.  2  of  Pig.  5b  and  could  be  substantia]  1 
reduced  if  two  Doppler  Signal  Processors  were  used  to  bypass  the  record¬ 
ing  of  one  of  the  heterodyne  signals. 

There  are  no  data  points  between  2  and  5  m  in  Pig.  S.  This 
region  coincides  with  the  maximum  of  the  scintillation  curve,  the  well- 
known  super-saturation  phenomenon.  Because  of  the  resulting  severe 
irradiance  fluctuations  that  prevail  under  these  conditions,  we  were 
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unable  to  make  the  proper  recording  and  processing  adjustments  that 
normally  minimize  and  stabilize  the  noise  level.  Several  measurements 
were  nevertheless  made  at  z  =  2.5  m,  but  the  results  were  unreproducible . 

In  spite  of  the  difficulties  created  by  the  relatively  large 

processing  noise,  Fig.  8  confirms  that  the  linear  relationship  given  by 

eq.  17  is  consistent.  The  slope  of  the  regression  line  can  now  be  used 

to  determine  the  value  of  the  average  velocity  U  of  the  convection  jets 

or  plumes.  U  is  the  important  but  yet  unspecified  parameter  of  eq.  15. 

The  measured  slope  is  4.0  x  10  ^  V2/ m  and  from  the  values  of  the  other 

parameters,  i.e.  K  =  10  V/rad/s,  k  =  9.93  x  106  rad/m, 

_4  -1/3 

C  =  1.05  x  10  m  ,  t  =  2  mm,  and  A  =  15  mm,  it  follows  from  eq.  17 
n  'o’  ’  M 

that 


U  =  27  mm/s. 


[18] 


This  value  is  in  excellent  agreement  with  visual  observations  of  the 
rising  convection  currents  in  the  turbulence  simulator.  Therefore, 
this  result  constitutes  further  evidence  that  eq.  15  is  a  justified 
approximation  or,  more  specifically,  that  the  measured  time  derivative 
of  the  turbulent  phase  is  indeed  linearly  related  to  the  vertical 
component  v  of  the  angle  of  arrival. 


UNCLASSIFIED 

26 


The  principal  objective  of  this  experiment  is  to  verify  the 
assumption  of  weak  correlation  between  the  angle  of  arrival  and  the 
irradiance.  Now  that  the  measured  d^/dt  has  been  proven  proportional 
to  the  vertical  component  of  the  angle  of  arrival,  we  can  use  the 
simultaneously  recorded  irradiance  and  phase-derivative  signals  to  test 
directly  this  hypothesis .  The  results  are  listed  in  Table  1. 


The  measurements  were  performed  with  the  Correlation  and  Probabi¬ 


lity  Analyzer.  The  cross -correlation  <v(I-<l>)>  was  estimated  by 

averaging  over  several  measurements  of  <(v  -v  )fl  -<1  >)>  and 

>'1  Xc  1  l 

<(v  -Vy  )(I  -<1  >)>  where  the  subscripts  1  and  2  refer  to  the  position 
points  r  and  r  in  the  transverse  plane.  We  find  that,  within  the 
experimental  errors,  both  expressions  give  approximately  the  same  value. 


thus  indicating  that  £  -  j:  is  sufficiently  Large  for  the  covariances 

<v  (I  -  <1  >)>  and  <v  (,1  -<I  >l>  to  vanish.  Finallv,  <\  •  was 
>'1  2  2  X;  ii'  X 

estimated  from  the  measured  using  the  theoretical  eq.  I" 

which  states  that 


-  |l  -  U  ‘8.1b  A‘/t‘]'1/b  }  _1 


-vv  )->. 

y 


where  A  =  jr  -  r  ,  =  15  mm.  The  experimental  error  is  relatively 
large  because  the  cross-correlation  value  is  almost  within  the  instru¬ 
mental  precision. 
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TABLE  I 


Comparison  of  the  measured  and  the  predicted  cross-correlations 
between  the  angle  of  arrival  and  the  irradiance. 


V  (I  - 

y 

<I>)>/^f  <v^><(I  -  <I>)Z>.  z: 

propagation  distance. 

i 

Z 

Y 

Y 

(m) 

(measured) 

(predicted,  eq.  22) 

1.20 

0.015  ±  0.01 

<  0.023 

1.62 

0.035  ±  0.04 

<  0.034 

3.00 

0.030  ±  0.04 

<  0.074 

3.47 

0.055  +  0.04 

<  0.087 

3.91 

0.070  +  0.04 

<  0.100 

4.40 

0.085  ±  0.04 

<  0.112 

The  results  of  Table  I  confirm  that  the  cross-correlation  between 
the  angle  of  arrival  ,v  and  the  irradiance  I  is  indeed  small,  actually 
smaller  than  8.5%  for  all  data  recorded  in  this  experiment.  Therefore, 
the  hypothesis  of  weak  correlation  between  v,  and  I,  which  is  basic 
to  our  theoretical  model  of  propagation  in  turbulence  (Refs.  1-5),  is 
validated  by  direct  measurements. 

From  the  theoretical  model  of  Refs.  2-3,  we  obtain  that 
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<v  (I-<I>)> 

<v2>1/2<fc-<l>)2>1/2 

y 

If  w(t)  is  the  radius  of  a  Gaussian  beam  in  turbulence,  and  if  the 
expression  of  Ref.  3  for  w(t)  is  used,  we  have 

<v  (I-<I>)>  -0 . 3S (C  /n  ) (kt  a<I>/o(v/wi 

<v2>1/2<(I-<I>)2>1//2  {w2+0.909(C  /n  )2k1/6;1  /b}1/2  <  ( I -<  I  >  J  ••>1 

y  o  no  v  ' 

[-1J 

where  wq  is  the  e-folding  radius  of  the  Gaussian  source  irradiance 
profile.  For  the  data  of  Table  I,  the  exact  coordinates  of  the  measure¬ 
ment  points  relative  to  the  beam  axis  were  unfortunately  not  recorded. 
Hence,  the  quantity  , 3<I>/9 (v/w) j /< (I -<l>) 2>  cannot  be  determined 

exactly  but  we  can  use  its  maximum  value,  equal  to  0.72  for  a  Gaussian 
,  1/2 

beam  with  <(I-<I>)‘;>  /<I>  =  1.2,  to  set  a  theoretical  limit.  It  thus 

follows  for  the  conditions  of  the  present  experiment  that 


<v  fl -<!>)> 1 


<  0.018  z" 


<V2,1/2<(I-<I>)2> ‘/:  I>  *  °-°»24  Z171*)1'2 

y  J 


where  the  propagation  distance  z  is  expressed  in  meters. 
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The  numbers  derived  from  eq.  22  are  compared  with  the  experimental 
data  in  Table  I.  The  ratio  of  the  measurements  to  the  theoretical  limits 
is  65%  on  the  average  which  constitutes  an  excellent  agreement  consider¬ 
ing  the  various  approximations  made  to  derive  the  simple  algebraic  ex¬ 
pression  of  eq.  22.  Therefore,  not  only  is  the  hypothesis  of  weak 
correlation  directly  verified  but  its  application  in  Refs. 2-3  to  the 
development  of  theoretical  expressions  for  the  cross-correlation  moments 
is  also  confirmed  by  direct  comparison  with  the  measurements  of  <v^I>. 
This  conclusion  constitutes  the  principal  result  and  fulfills  the  main 
objective  of  the  present  study. 

3.3  Probability  Distribution  of  the  Angle  of  Arrival 

Another  important  hypothesis  made  in  the  development  of  our  prop¬ 
agation  model  concerns  the  probability  distribution  of  the  angle  of 
arrival  v.  It  was  assumed  that  ;v  obeys  a  normal  or  Gaussian  probability- 
density  function  which  was  used  to  relate  the  third-  and  the  fourth- 
order  statistical  moments  of  to  its  second-order  moments.  This  hypo¬ 
thesis  can  now  be  tested  directly  from  the  recorded  data  of  the  present 
experiment . 

The  probability  density  functions  of  the  difference  between  the 

demodulated  heterodyne  signals,  i.e.  (e(r  )  -  e(r  )),  were  calculated 

1  2 

at  different  propagation  distances  with  the  Correlation  and  Probability 
Analyzer.  The  instrument  provides  a  100-point  analysis  of  the  density 
function.  For  all  cases  reported  here,  the  electrical  voltage 
(e(r^)  -  e^r^))  was  sampled  at  a  rate  of  500  Hz,  and  2^  consecutive 
points  were  used.  A  typical  measured  probability  function  is  plotted 
in  Fig.  9a.  Comparison  with  a  least-square  fitted  Gaussian  curve  shows 
that  the  measured  function  is  not  Gaussian  as  illustrated  in  Fig.  9a. 


Probability  Probability 
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We  have  already  discussed  that  the  measured  random  voltage 
difference  (e(r  )  -  e(r  ))  comprises  a  significant  noise  contribution. 

To  evaluate  the  latter,  the  frequency-demodulated  signal  from  one  detec¬ 
tor  was  recorded  simultaneously  with  the  heterodyne  signal  from  the  same 
detector.  On  playing  back  the  tape,  the  difference  (e(r^)-e(j:  ))  between 
the  recorded  demodulated  output  and  the  demodulated  output  of  the  re¬ 
corded  heterodyne  signal,  which  constitutes  a  sample  of  the  noise  contri¬ 
bution  to  (e(x  )-e(£  )),  can  be  obtained.  For  example,  the  probability 
density  function  of  the  noise  sample  corresponding  to  the  data  of  Fig. 

9a  is  drawn  in  Fig.  9b.  It  is  evident  from  the  observation  of  these  two 
figures  that  the  standard  deviation  of  the  noise  is  not  negligible 
compared  with  that  of  the  combined  signal -plus -noise.  Therefore,  a 
deconvolution  is  necessary  to  recover  the  true  probability  function  of 
the  phase-derivative  signal.  Indeed,  since  the  noise  and  the  signal  are 
uncorrelated,  the  measured  probability  function,  pe(x),  is  given  by  the 
convolution  of  the  respective  probability  functions  of  the  signal, 
ps(x),  and  of  the  noise,  pb(x),  i.e. 


pew 


ps(*-y)  pb(y)  dy  , 


[23] 


where  pfi(x)  is  defined  as  the  probability  that  the  measured  voltage 
(e(r  )  -  e(x  ))  lies  between  x  and  x  +  dx,  and  similarly  for  p  (x) 

12  S 

and  Pb(x) . 

The  required  function  pg (x)  is  derived  by  an  iterative  method. 
The  approach  consists  in  setting  pg(x)  =  A(x)  exp(-yx2)  and  determining 
A(x)  and  y  by  successive  approximations  of  the  convolution  integral  of 
eq.  23.  The  method  was  found  stable  for  all  cases  analyzed  in  this 
experiment  and  it  converges  after  a  few  iterations.  The  demodulated 
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probability  function  pg(x)  from  the  experimental  data  of  Figs.  9a-b  is 
plotted  in  Fig.  9c.  By  substituting  this  function  pg(x)  back  into 
eq.  23,  we  can  verify  the  accuracy  of  the  method.  This  result  is 
illustrated  in  Fig.  9d  where  the  calculated  and  the  experimental  pg(x) 
are  compared.  The  agreement  is  excellent,  the  root -mean -square 
difference  is  less  than  2%  of  the  maximum  probability,  which  is  of  the 
order  of  the  experimental  scatter. 

The  most  straightforward  test  for  Gaussian  or  normal  statistics 
consists  in  plotting  the  experimental  probability  distributions  on  a 
normal  probability  scale.  This  is  done  in  Fig.  10  for  all  the  measured 
deconvoluted  probability  distributions.  If  the  normal  hypothesis  is 
correct,  the  data  must  fall  on  a  straight  line,  which  is  well  verified 
in  Fig.  10  over  the  range  spanning  from  -4o  to  4o,  where  o  is  the  standard 
deviation.  It  is  worth  emphasizing  that  the  interval  ±  4o  encloses  the 
random  signal  with  a  probability  of  99.98°o.  Therefore,  the  hypothesis 
of  normal  statistics  for  (vy(r^)  -  v^Qr  ))  is  fully  justified  on  the 
basis  of  this  test. 


A  quantitative  measure  of  the  degree  of  deviation  from  the 
hypothetical  normal  distribution  is  provided  by  the  well-known  chi-square 
test  (Ref.  12()  .  The  test  is  based  on  the  fact  that  for  samples  large 
enough  the  quantity 


(f.  -  Mp.)2 
Mp. 


1 24] 


has  approximately  the  chi-square  distribution  with  e  =  m  -  1  -  s 
degrees  of  freedom,  f.  denotes  the  number  of  counts  or  events  registered 
in  the  it^1  channel,  p.  is  the  tested  normal  probability  function,  i.e. 


Probab i I t  ty  <  %  ) 
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FIGURE  10  -  Deconvoluted  probability  distributions  of  the  difference  in 
the  angle  of  arrival  plotted  on  the  normal  probability  scale 
versus  the  centered  signal  amplitude  normalized  by  the 
standard  deviation  o.  The  normal  distribution  is  represented 
by  the  straight  line.  7:  z  =  1.20  m;  A:  z  *  1.62  m; 

O:  z  =  3.00  m;  0:  z  =  3.47  m;  O:  z  =  3.91  m;  x:  z  =  4.40m. 


bn 
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Pi 


dx  exp  [ - (x-u) 2/2c2 ] 


[25] 


M  is  the  size  of  the  sample,  and  s  is  the  number  of  parameters  estimated 
from  the  sample  which  in  this  case  is  equal  to  2,  the  expectation 
value  u  and  the  standard  deviation  a.  The  test  consists  in  calculating 
X2  according  to  eq.  24  and  comparing  the  values  to  the  percentiles  of 
the  chi-square  distribution. 

The  probability  data  were  calculated  from  2^  consecutive  points 
sampled  at  a  rate  of  500  Hz.  Since  the  correlation  time  of  the  signal 
is  typically  equal  to  30  ms,  the  number  of  independent  sampled  events  is 
about  M  —  2  /15  =  8738.  The  number  of  channels  provided  by  the  instru¬ 

ment  is  100  but  regrouping  was  performed  to  satisfy  the  generally  accepted 
condition  f^  >  10  for  the  test  to  be  significant.  The  results  are  listed 
in  Table  II.  For  all  cases,  we  find  that  x2  <  X297  where  X2t>7  defined 
such  that 


P(X2  >X297)  =  97*, 


[26] 


where  P  is  the  probability  calculated  from  the  chi-square  distribution. 
Therefore,  the  test  shows  that  the  normal  probability  hypothesis  is 
acceptable  on  the  97°0  level,  i.e.  with  a  97“  probability  of  rejecting 
the  hypothesis  when  in  fact  it  is  true,  which  constitutes  a  very 
stringent  test. 
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In  summary,  the  model  of  normal  distribution  for  the  statistics 

(v  (r_  )  -  v  (r))  is  valid  with  a  high  degree  of  accuracy.  Since 

Vy(r^)  and  v  Cl  )  are  nearly  independent  because  of  the  large  relative 

separation  distance  ir  -  x  |  =  IS  mm,  and  since  the  convolution  of  two 

~] 1  2 

Gaussian  functions  gives  a  Gaussian  function,  it  is  highly  plausible 
that  the  individual  v  (r  )  and  v  (r  )  are  also  normally  distributed. 


A 


TABLE  II 


Results  of  the  chi-square  test  on  the  normal  probability  distribution 
hypothesis  for  the  turbulent  angle  of  arrival.  x2»  chi-square  parameter 
calculed  according  to  eq.  24;  X2gy»  97%  -  percentile  of  the  chi-square 
distribution;  M,  sample  size;  l,  number  of  degrees  of  freedom;  z, 
propagation  distance. 
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Moreover,  because  no  asymmetry  has  been  observed  in  the  statistical 
properties  of  the  laser  beam,  it  is  equally  probable  that  the  horizontal 
component  v  (r)  is  also  normally  distributed.  Hence,  we  conclude  that 
the  hypothesis  of  normal  probability  distribution  of  the  angle  of  arrival 
v,  used  in  the  development  of  our  propagation  model,  is  experimentally 
confirmed  with  a  high  degree  of  confidence. 


4.0  DISCUSSION  AND  CONCLUSION 


The  application  of  the  heterodyne  method  to  measure  the  instan¬ 
taneous  turbulent  fluctuations  of  the  angle  of  arrival  led  to  some  exper¬ 
imental  difficulties.  There  was  the  problem  of  sensitivity  to  mechanical 
vibrations  which  was  circumvented  by  simultaneous  detection  at  two  points 
and  subsequent  subtraction  of  the  recorded  signals.  This  made  it  much 
more  difficult  to  achieve  and  maintain  the  proper  alignment  of  the  inter¬ 
fering  beams.  However,  the  main  source  of  error  came  from  the  very 
nature  of  the  detected  phenomenon.  Indeed,  the  heterodyne  signals  are 
amplitude  modulated  by  the  turbulence  over  a  dynamic  range  wider  than 
that  of  the  magnetic  tape  recorder.  Consequently,  the  recorded  hetero¬ 
dyne  signal  frequently  falls  below  the  noise  level  _f  the  recorder  which 
causes  random  erroneous  peaks  in  the  output  of  the  frequency  demodulator. 
Coupled  with  the  relatively  small  voltage-frequency  sensitivity  of  the 
instrument  used,  this  produced  a  low  signal -to-noise  ratio.  Hence, 
careful  attention  had  to  be  devoted  to  all  experimental  and  analytical 
steps  to  extract  the  useful  information  from  the  output  signal  of  the 
frequency  demodulator.  Meaningful  results  were  nevertheless  obtained 
since  the  noise  and  phase  contributions  are  mutually  uncorrelated. 
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Despite  the  difficulties  encountered,  our  experiment  was  very 
successful.  It  demonstrated  that  the  time  derivative  of  the  phase 
measured  from  the  frequency-demodulated  heterodyne  signal  is  indeed 
linearly  related  to  the  vertical  component  of  the  angle  of  arrival. 

This  permitted  the  calculation  of  the  statistical  correlation  between  the 
angle  of  arrival  and  the  irradiance.  The  hypothesis  of  weak  correlation, 
which  is  central  to  the  development  of  our  model  of  propagation  in 
turbulence  (Refs.  1-3),  was  thus  directly  confirmed.  Also  verified 
were  the  resulting  theoretical  predictions  for  the  cross-correlation 
moments.  Finally,  the  assumption  of  normal  probability  distribution 
for  the  angle  of  arrival,  which  was  also  used  in  the  construction  of 
our  propagation  model,  is  in  excellent  agreement  with  the  direct  proba¬ 
bility  data  of  this  experiment. 

In  conclusion,  this  report  confirms,  by  comparison  with  measure¬ 
ments,  two  fundamental  assumptions  of  the  DREV  propagation  model  in 
turbulence.  Of  course,  these  had  already  been  proven  consistent  since 
the  theoretical  predictions  of  Refs.  1-3  regarding  the  average  ir-adiance 
and  the  irradiance  variance  profiles  of  laser  beams  traveling  in 
turbulent  media  were  verified  over  a  wide  range  of  turbulence  and  beam 
conditions.  The  merit  of  the  present  experiment  consists  of  the 
direct  verification  of  these  two  hypotheses  for  conditions  ranging  from 
the  weak  to  the  strong  scintillation  limits. 
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